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We show that the strong hybridization between the iron 3d and the arsenic 4p orbitals, in the newly 
discovered iron-based high-T c superconductors, leads to an explanation of certain experimental 
observations that are presently not well understood. The existence of a lattice distortion, the 
smallness of the Fe magnetic moment in the undoped systems, and the suppression of both the 
lattice distortion and the magnetic order upon doping with fluorine, are all shown to result from 
this hybridization. 



An interetsing development occurred recently with the 
discovery of a new class of layered, iron-based, high tem- 
perature superconductors. Kamihara et alP reported a 
superconducting transition temperature T c =26 K in fluo- 
rine doped LaOFeAs, and T c increased under pressure to 
43 K. 2 Replacement of lanthanum with other rare earth 
elements gave a series of superconducting compounds 
ReOi-^F^FeAs where Re = Ce, Pr, Nd, or Sm , with 
transition temperatures close to or exceeding 50 KJ 3 * 4 * 5 * 6 *. 
Oxygen-deficient samples were also fabricated and found 
to superconduct at 55 K.™ The parent compound, Re- 
OFeAs, has a tetragonal unit cell with two ReOFeAs 
molecules, and it consists of a stack of alternating ReO 
and FeAs layers. The FeAs layer, in reality, consists of a 
square planar sheet of Fe atoms sandwiched between two 
sheets of As atoms, such that each Fe atom is tetrahe- 
drally coordinated with four As atoms. 

Certain facts that emerged from neutron scattering 
measurements are not well understood: 

1) A tetragonal-to-orthorhobic structural phase tran- 
sition takes place in the undoped system at Ts^l50 K, 
below which a n antifcrromagnetic (AFM) order is estab- 
lished mmm 

It was suggested^ that the tetragonal-to- 
orthorhombic distortion results from magnetic frustra- 
tion. However, the magnetic order oc curs more than 20 
K below T s in LaOFeAs and PrOFeASj^casting doubt 
on such an interpretation. 

2) The me asured m agnetic moment of Fe is very small, 
~0.4-0.8 >u, ^ j S | 9 | io |TT] wnereas calculations based on den- 
sity functional theory (DFT) give -2 ^GUMMED 

3) It is observed that upon replacing a small fraction of 
the oxygen atoms in the ReO layer by fluorine atoms, the 
tetragonal-to-orthorhombic distorti on is gr adually sup- 
pressed, along with the AFM order PEESffi] 

In this letter we present a scheme within which the 
above observations may be explained. It is based on the 
consequences of the strong hybridization between the Fe 
3d and the As 4p orbitals. On the one hand, this hy- 
bridization gives rise to superexchange between different 
Fe atoms, leading to the establishment of the AFM or- 
der. On the other hand, the hybridization will have a 
strong effect on the states of any given Fe atom; it is this 
effect that we explore here, and show that it offers an 
explanation of the above observations. 
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FIG. 1: The occupation scheme of the Fe 3d levels in the 
a) tetragonal and b) orthorhombic structures of ReOFeAs, 
where Re is a lanthanide. 

In the tetragonal phase, and in the absence of hy- 
bridization, the crystal field at a given Fe site, pro- 
duced by the four surrounding As ions, has S4 symme- 
try, and it splits the otherwise degenerate Fe 3d levels 
into three nondegenerate levels and one doubly degener- 
ate levefS3 as shown in Fig. [lp,. Since DFT calculations 
indicate that an Fe site has a total spin 1, the six 3d 
electrons are taken to occupy orbitals as shown in Fig. [T| 
this occupation scheme was also suggested by LP^ and 
Baskaran 18 . Following Li,E3 we assume that the energy 
splitting produced by the crystal field is sufficiently large 
to neglect the filled 3d 2 2 and 3d xy orbitals and the empty 
3d x 2_ y 2 orbital in discussing the low energy behavior. A 
tetragonal-to-orthorhombic lattice distortion will lift the 
remaining degeneracy, and assuming that the lattice con- 
stant a is larger than b, the energy level scheme of the d 
orbitals will be as shown in Fig. [Tb. If in Fig. [T^, we set 
e d^ = td yz = e, then in Fig. lib we have e diz = e-A = ei 
and €d yz — e + A = £2- We consider a system compris- 
ing an Fe atom and an adjacent As atom, for which the 
Hamiltonian may be expressed as 

H = Ti-o + f~i , "Ho = 7~(-d + Hp, 
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Using H = — [Hq , S] . we obtain 



FIG. 2: The state |1) along with the other four states that 
can be reached from it through hybridization of the Fe 3d and 
the As 4p orbitals. 



Here Tid and H p are the Hamiltonians for the Fe and 
As atoms, respectively, c pr7 (c pcr ) is the creation (anni- 
hilation) operator for an electron in the As 4p orbital 
with spin a, c' d a {cd xzi a) creates (annihilates) an elec- 
tron with spin a in the Fe 3d xz orbital, with a simillar 
definition applying to c d .(cd BZ , CT ), a and (3 are hopping 
matrix elements, taken for simplicity to be constant and 
real, and H.C. stands for hermitian conjugates of the last 
two terms in Ti . The terms in Ti express the hybridiza- 
tion between the Fe 3d orbitals and the As 4p orbitals, 
in the form of electron hopping between these orbitals. 

The states shown in Fig. [l]are eigenstates of Ho, but 
not H. Nonvanishing matrix elements connect the state 
|1) to other states shown in Fig. |j State | 2) can be 
reached from |1) by the spin- up electron in Fe 3d yz hop- 
ping to the As ion, followed by a spin-down electron hop- 
ping from the As ion to the Fe 3d xz orbital. Alternatively, 
a spin-down electron may hop first from the As ion to the 
Fe 3d xz orbital, followed by the electron in the Fe 3d yz 
orbital hopping to the As ion. The ground state of the 
Fe atom is now obtaind by diagnolizing the Hamiltonian 
matrix in the subspace spanned by these five states. 

The Hamiltonian is projected onto this 5-dimensional 
subspace by a canonical transformation 
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Choosing the operator S such that [H ,S] = — H , and 
keeping terms to 2 nd order in H , we obtain 

Ht — Ho + H e ff 
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The eigenstates |m) of Ho are in fact composite states 
comprising the 3d xz and 3d yz orbitals of the Fe atom, 
along with the 4p orbitals of As. We carried out DFT 
calculations on undoped LaOFeAs, along with Bader's 
atoms in molecules (AIM) theorypa and found that As 
has a charge of — 0.92e and Fe has a charge of 0.30e, where 
e is the charge of the proton. The charge transferred from 
LaO to the FeAs layer resides mainly on the As 4p 2 or- 
bital, which has an energy much lower than that of the 
As 4p x or 4p y orbital. Furthermore, by examining the 
projected density of states (PDOS) reported by Cao et 
alP^we see that the DOS at the Fermi level is dominated 
by the Fe 3d xz and 3d yz states, while the DOS of the p 
states has 3 peaks: one peak, which results from filled 
4p z bands, is 3-4 eV below the Fermi energy E*, and 2 
other peaks, one being about 0.5 eV below and the other 
1 eV above Ep; these two peaks result from the 4pa; and 
4p y orbitals. Since there are two As atoms per unit cell, 
we can assume that there are two occupied 4p z bands, 
3-4 eV below E^ , two occupied bands with mixed 4p.,; 
and 4p y character, ~0.5 eV below Ejr, and another two 
empty bands with mixed 4p2, and 4p a character, ~1 eV 
above E^. Considering any single As ion, we may assume 
that it has four electrons in the 4p orbitals, two of which 
occupy the 4p 2 orbital, and the other two electrons oc- 
cupy an orbital with mixed 4p;,; and 4p a character which, 
for convenience, we simply call p, and there is an empty 
orbital, also with mixed 4jp x and 4pj, character which we 
call p'. We ignore the filled, low energy 4p z orbital , and 
we describe the state of the As atom by the occupation of 
the p and p orbitals. The eigenstate |1) of Ho is written 
explicitly as 

|1) = \Fe;d xzA ,d yzA )\As;p u ). 

The eigenstate |2) of Ho is 

|2) - \Fe;d xzAl )\As;p h p\) 

Note that |2) differs from |1) by having both d electrons 
in the d xz orbital and by a spin flip in the As p orbitals. 
In order to construct the matrix element (2|iJ e *Jl) we 
need the intermediate states \k). As discussed earliear 
there are two such states, 

= \Fe;d xzA }\As; P n,p |), 



where Heff — \\H ,S]. The matrix elements of H e jy, 
connecting eigenstates of Hq, are given by 



(n\H eff \m) = -(n\HS-SH\m) 



= lY^ n \H'\k)(k\S\m) - (n\S\k)(k\H'\ 



\k 2 ) = \Fe;d xz H,d yz T)|As;pT) 

Because of the onsite Coulomb repulsion on the d orbital, 
the state \k 2 ) has s much higher energy than \ki), and it 
will be ignored. We can then write 
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where J is the Hund's exchange term in the d orbitals, 
e p (e p ') is the energy of an electron in a p(p ) orbital on 
the As atom, and Hund's exchange term in the As 4p 
orbitals is ignored. The factor of 4 results from the fact 
that an electron on the Fe atom can hop to any of the 
four surrounding As atoms. In a similar way, we can 
construct the remaining matrix elements and diagonalize 
the Hamiltonian in the 5-dimensional subspace spanned 
by the states |1), |5). 

Needless to say, the expressions for the eigenvalues and 
eigenvectors of the 5x5 matrix are unwieldy. Instead, 
we diagonalize the matrix in the 2-dimcnsional subspace 
spanned by |1) and |2), draw the necessary conclusions, 
and then extrapolate to the more general case. 

Thus, assuming for the time being that only states |1) 
and 1 2) are mixed by the Hamiltonian, we can express 
Ht by a 2 x 2 matrix 



ei + e 2 - J + 2e p 7 

7 2ei + e p + e p > 

The ground state energy is now given by 



1, 

E gs = ei + e p + -{ex + e 2 + e p > + e p - J) 



\^(ei-e 2 + e p , - e p + J) 2 + 4| 7 | 2 



and E gs is lower than the energy E\ of state |1) by 5E 
given by 

SE = E gs -E 1 = i[(J-2A + 2A p ) 



where 2A 



p 



(J-2A + 2A P ) 2 +4| 7 | 2 ] 



e p . The unperturbed energy Ei 



ei + 62 — J + 2e p — 2e + 2e p — J is the same in both the 
tetragonal and orthorhombic phases, but 6E does depend 
on the crystalline phase. Since 7 depends only weakly on 
A, we have 
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where the superscript e stands for 'electronic'. The above 
expression is negative for all values of J, A and 7. From 
our earlier discussion, we may take — e p ~ 0.5 eV, 
( d — e p ' ~ —1 eV, 2A p ~ 1.5 eV. The parameters J, 
a, and j3 are all estimated to be in the range ~0.7-l 
eV* 16 * 17 ^ this gives 7 ~ 1.5 — 3 eV. The above expression 
then reduces to 



5E { ;X ho -5E tetra 



(e) _ 



-A(l- 
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which is negative. We thus conclude that a tetragonal- 
to-orthorhombic distortion lowers the electronic energy 



of the d state by a term proportional to A (ignoring the 
A 2 term, being smaller than the linear term). The dis- 
tortion takes place by a small increase in the lattice con- 
stant a and a corresponding decrease in lattice constant 
b. Assuming this small distortion is Sx , then A is pro- 
portional to Sx, and we can write SE^ = £,5x, where £ is 
some constant. On the other hand, this distortion raises 
the elastic energy by an amount 5E L = rj(5x) 2 where r\ 
is some constant. The value of the distortion Sx is found 
by minimizing 



SE 



total 



5E^ e 



SE 
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giving Sx — £/2?y, and SE total = — ^ 2 /4rj. What we have 
just shown is that the hybridization between the Fe 3d 
states and As 4p states causes a greater downward shift 
in energy if the crystal distorts from tetragonal to or- 
thorhombic, and seeking the lowest energy available, the 
crystal will indeed undergo such a distortion. 

The ground state eigenket is obtained from the eigen- 
vector corresponding to the ground state energy, 



|Vv) = Ci|l) + c 2 |2), 



where 
ci 



2| 7 



2A P - 2A) 2 + 4 7 2 - ( J + 2A P - 2A)] . 



Using the values of the parameters given earlier, we find 
\cijc\\ in the range ~0.5-0.7. In the ground state the 
expectation value of the magnetic moment of the Fe atom 
is 

( Ai ) = | Cl | 2 ( M ) 1 + |c 2 | 2 ( M ) 2 

Since (/x)2, the magnetic moment in state |2), is zero, 
it follows that (/i) ~ 0.7(/i)i. In other words, the hy- 
bridization causes a large reduction in the value of the 
would-be magnetic moment. We may now consider the 
mixing of all the states shown in Fig. [2] The qualitative 
results obtained so far remain unchanged. Qualitatively, 
the ground state will be a mixture of five states, in four 
of which the magnetic moment is zero. Therefore, we ex- 
pect a further reduction in the expectation value of the Fe 
magnetic moment. In fact we do not need to restrict our- 
selves to these states; we may also consider mixing with 
states where one electron resides on the A x i_ y i orbital. 
All this leads to the conclusion that strong hybridization 
between the Fe 3d and As 4p orbitals causes spin fluctu- 
ations on the Fe sites and a subsequent reduction in the 
measured value of the Fe magnetic moment , in agreement 
with neutron scattering measurements.^ 4 * 5 * 6 ! 

We now consider what happens as some oxygen atoms 
are replaced by fluorine atoms. Fluorine doping causes 
further negative charge transfer from the ReO /F layer to 
the As atoms, since it involves replacing 2 ~ with F~. 
This strengthens the ionic bond and shortens the distance 
between the Re and As ions. Indeed, neutron scattering 
measurements El indicate that the La- As separation de- 
creases by 0.03 A upon F-doping, with a concomitant 
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reduction in the c-lattice constant. But the extra charge, 
donated by the Re ions, is not uniformly spread among 
the As ions; more of this charge resides on As ions close 
to the F atoms, rather than on those further away, and 
some local distortion of the As ion positions is expected. 
With such distortion, and with an uneven charge distri- 
bution on the As ions, the crystal field at a given Fe site 
now has a symmetry lower than S±. Considering a single 
Fe atom, the Hamiltonian is written as 

H = Ho + V 

where H.q contains the part of the crystal field with S4 
symmetry, while V is the additional part which lowers 
that symmetry. The orbitals Fe 3d xz and 3d yz are eigen- 
states of H.Q but not TL; the perturbation V mixes these 
two states. Let (3d xz \V\3d xz ) — Vn, (3d yz \V\3d yz ) = 
V22, and (3d xz \V\3d yz ) = V12; these matrix elements are 
proportional to the F-dopant concentration. Then the 
energy splitting 2A of the two states is given by 

2A= yVn - v 22 ) 2 + 4Vi 2 2 , 

and the two new split states \a) and |&) arc linear com- 
bination of \3d xz ) and \3d yz ). As the F dopant concen- 
tration increases, the energy splitting increases, and if 
2A becomes larger than J, the contribution of state |2), 
shown in Fig [2] to the ground state wave function be- 
comes dominant, and the Fe magnetic moment becomes 
very small, causing the magnetic order to be suppressed. 
If the fluorine concentration is x, then in the ReO/F 
layer, every position in the O/F plane is occupied by 
a fluorine atom with probability x and by an oxygen 
atom with probability 1 — x. The overall symmetry of 
the crystal structure, as seen in X-ray measurements, is 
still tetragonal with symmetry group P4/nmm. The in- 
stantaneous potential seen by an electron depends on the 
ionic positions and is given by 

V = V s + Y / (dV/dQ k )Q k + ... 

k 



where V s is the static part of the potential, Qfc is a lattice 
normal coordinate, and because of the invariance of V 
under the crystal symmetry group, dV/dQk transforms 
according to the same irreducible representation as 
Since in the presence of F doping, the Fe electronic or- 
bitals \i) are nondegenerate, then (i\dV / dQk\i) vanishes 
unless dV/dQk is totally symmetric under the crystal 
symmetry group operations. It follows that any sponta- 
neous distortion that occurs must be totally symmetric 
and will not lower the crystal symmetry. Hence, F doping 
suppresses the tetragonal-to-orthorhombic distortion. 



In conclusion, we have shown that strong hybridization 
between the Fe 3d orbitals and the As 4p orbitals intro- 
duces some interesting effects. In the undoped systems, it 
causes a larger downward shift in the electronic energy if 
the crystal undergoes a tetragonal-to-orthorhombic dis- 
tortion which lifts the degeneracy of the Fe 3d xz and 3d yz 
levels. The ground state, in the presence of hybridization, 
is found to be a mixture of states, most of them with zero 
Fe magnetic moment, so that in this ground state the ex- 
pectation value of the Fe magnetic moment is small. We 
have also argued that the symmetry lowering caused by 
fluorine doping suppresses both the lattice distortion and 
the magnetic order of the Fe moments. Finally, we note 
that although the mechanism for superconductivity in 
the iron-based superconductors is not understood yet, we 
believe that the spin fluctuations on the Fe sites, caused 
by hybridization between the Fe 3d and the As 4p states, 
may play a role, not only in suppressing the magnetic 
order, but also in the emergence of superconductivity in 
these compounds. 
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